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I.  Executive  summary  of  accomplishments  relative  to  SOW 


In  the  Statement  of  Work  (SOW)  mutually  agreed  upon  with  DARPA,  we  were  tasked 
with  the  following: 

•  Task  1.  Search  for  new  quantum  algorithms 

•  Task  2.  Improve  theoretical  understanding  of  existing  quantum  algorithms 

•  Task  3.  Analyze  issues  associated  with  implementation 

In  regard  to  Task  1,  we  have  accomplished  the  following: 

•  We  have  created  six  new  quantum  hidden  subgroup  (QHS)  algorithms: 

o  A  continuous  Shor  algorithm 
o  The  wandering  Shor  algorithm 
o  The  lifted  Shor  algorithm 
o  A  quantum  circle  algorithm 
o  A  QHS  algorithm  dual  to  Shor’s  algorithm 
o  A  QHS  algorithm  for  functional  integrals 

In  regard  to  Task  2,  we  have  accomplished  the  following: 

•  In  creating  the  above  algorithms,  we  have  found  three  new  systematic  procedures 
for  creating  and  finding  new  quantum  algorithms,  namely: 

o  Lifting 
o  Pushing 
o  Duality 

•  By  applying  the  methods  of  combinatorial  group  theory,  we  have  made 
substantial  progress  in  the  theoretical  understanding  and  analysis  of  non-abelian 
quantum  hidden  subgroup  (QHS)  algorithms.  In  particular, 

o  We  have  shown  that  non-abelian  QHS  algorithms  can  be  lifted  to  quantum 
algorithms  on  free  groups  in  polytime,  and  then  pushed  back  to  their 
original  domain  also  in  polytime.  Hence,  there  is  no  loss  of  generality 
resulting  from  focusing  only  on  non-abelian  QHS  algorithms  on  free 
groups,  where  the  technical  problems  are  substantially  simpler.  This  is  a 
significant  result. 

o  We  also  have  shown  that  Schreier  2-sided  transversals  provide  a  best  way 
to  lift  non-abelian  QHS  algorithms  to  free  groups.  This  is  a  generalization 
of  what  we  call  the  Shor  transversal,  which  was  used  to  create  the 
wandering  Shor  algorithm. 

•  We  have  demonstrated  that  Grover’s  and  Shor’s  quantum  algorithms  are  more 
closely  related  than  previously  thought. 
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o  We  have  shown  that  Grover’s  algorithm  is  a  QHS  algorithm  in  the  sense 
that  it  solves  a  non-abelian  quantum  hidden  subgroup  problem. 

o  We  have  also  shown  that  this  QHS  problem  can  not  be  solved  by  the 
standard  non-abelian  QHS  algorithm. 

•  We  have  used  the  methods  and  techniques  of  quantum  topology  to  obtain  new 
results  in  quantum  computing.  In  particular, 

o  We  have  found  all  4x4  unitary  solutions  of  the  Yang-Baxter  equation,  and 
determined  which  of  these  are  universal  quantum  computing  gates. 

o  We  have  found  relationships  between  quantum  entanglement  and 
topological  linking. 

o  We  have  determined  new  universal  quantum  gates  that  are  solutions  with 
spectral  parameters  of  the  Yang-Baxter  equation. 

o  We  have  shown  how  quantum  teleportation  can  be  understood  in  terms  of 
the  categorical  formalism  of  quantum  topology.  This  approach  shows 
promise  of  being  a  useful  tool  for  distributed  quantum  computing. 

o  We  have  found  a  new  approach  to  creating  unitary  representations  of  the 
braid  group  based  on  the  bracket  model  of  the  Jones  polynomial.  This 
approach  includes  Kitaev’s  Fibonnaci  model,  and  shows  promise  of 
yielding  new  insight  into  anyonic  topological  quantum  computation. 

o  We  have  classified  all  representations  of  the  three-strand  braid  group. 

o  We  are  currently  writing  a  paper  that  shows  that  the  quantum  algorithm 
created  by  Aharonov,  Freedman,  Kitaev,  Jones,  and  Landau  for  finding  the 
values  of  the  Jones  polynomial  can  not  successfully  be  extended  via 
polynomial  interpolation  to  a  quantum  algorithm  that  actually  computes 
the  full  Jones  polynomial. 

o  We  have  developed  a  new  combinatorial  criterion  for  determining  whether 
or  not  a  pure  state  is  entangled. 

In  regard  to  Task  3,  we  have  accomplished  the  following: 

•  We  have  proposed  distributed  quantum  computing  (DQC)  as  a  roadmap  to 
scalable  quantum  computing,  i.e.,  as  a  strategy  for  effectively  using  existing  or 
near  future  (within  five  years)  quantum  computing  devices  (hardware)  to  perform 
BIG  tasks,  normally  thought  only  possible  on  non-existent  large  quantum 
computers.  By  DQC,  we  mean  quantum  computing  on  a  network  of  small 
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quantum  computers  interconnected  by  quantum  Einstein  Poololsky  Rosen  (EPR)  and 
classical  channels.  Any  existing  quantum  computer  device  which,  for  example,  can  transform 
photon  (flying)  qubits  into  system  qubits  and  back  (such  as  ion  traps,  neutral  atom 
devices,  linear  optics,  etc.)  could  be  used  to  form  such  a  quantum  network.  We 
have  obtained  the  following  results  in  DQC: 

o  We  have  created  a  universal  set  of  quantum  distributive  computing 
primitives  for  transforming  existing  quantum  algorithms  into  distributed 
quantum  algorithms,  namely: 

■  Cat-Creator 

■  Disentangler 

■  Reset 

■  Swap-Reset 

o  Based  on  the  above  DQC  primitives,  we  have  created  a  systematic 
methodology  for  transforming  non-distributed  quantum  algorithms  into 
distributed  quantum  algorithms.  This  procedure  can  be  automated. 

o  We  have  used  the  above  mentioned  systematic  methodology  to  create  a 
distributed  version  of  the  quantum  Fourier  transform  and  of  Shor’s 
quantum  factoring  algorithm. 

o  We  have  shown  that  distributed  quantum  algorithms  are  not  significantly 
less  efficient  than  their  non-distributed  counterparts.  In  particular,  we 
have  shown  that  the  increase  in  computational  complexity  (resulting  from 
distributed  computing  overhead)  is  insignificant  in  comparison  with  the 
greater  efficiency  and  speed  of  quantum  algorithms. 

o  We  have  also  pointed  out  that  DQC  provides  a  mechanism  for  better 
dealing  with  the  problem  of  decoherence,  i.e.,  it  provides  a  “divide  and 
conquer”  strategy  for  dealing  with  decoherence.  Once  EPR  channels  have 
been  established,  one  need  only  focus  on  the  decoherence  problem  for 
each  spatially  separated  quantum  device  in  the  network.  The  key  idea  is 
that  NOT  all  environmentally  entangling  transformations  are  equally 
likely.  In  particular,  for  spatially  separated  physical  quantum  computing 
devices,  the  most  likely  environmentally  entangling  transformations  are 
those  that  are  isolated  to  the  local  quantum  device  and  its  immediate 
environment.  (This  observation  is  debated  by  some  in  the  quantum 
computing  community.) 

Also  in  regard  to  all  Tasks  1,  2, 3,  we  have  accomplished  the  following: 

•  Asa  research  tool,  we  have  developed  a  Maple  software  package  for  simulating 
quantum  systems.  This  software,  as  well  as  an  instruction  manual,  can  be  found 
at  the  website: 

http://www.csee.umbc.edU/~lomonaco/DARPA/Q  1  -06finalrpt 
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•  Other  research  related  to  all  tasks,  but  either  not  yet  completed  or  not  yet 
undertaken,  is  the  following: 

o  The  use  of  QHS  algorithms  to  extract  geometric  patterns  from  photos. 

o  The  design  of  DQC  algorithms  for  systolic  arrays  of  quantum  smidgets, 
where  by  a  quantum  smidget,  we  mean  a  small  but  powerful  quantum 
computing  device,  such  as  a  quantum  linear  feedback  shift  register. 

II.  Accomplishments 

II.A.  Contributions  to  quantum  algorithms 

II.A.l.  Definition  of  quantum  hidden  subgroup  algorithms 

What  is  a  hidden  subgroup  problem?  What  is  a  hidden  subgroup  algorithm? 

Definition.  A  map  (p:G  — >  S  from  a  group  G  into  a  set  S  is  said  to  have  hidden 
subgroup  structure  if  there  exists  a  subgroup  K  of  G ,  called  a  hidden  subgroup,  and 
an  injection  i  :  G  /  K  — »  S ,  called  a  hidden  injection,  such  that  the  diagram 

G  — ^  S 

v\  A 

G'Ky 

is  commutative1 ,  where  G  /  K  denotes  the  collection  of  right  cosets  of  K  in  G ,  and 
where  v  :G  — >  G! K  is  the  natural  surjection  of  G  onto  Gl K  .  We  refer  to  the  group 
G  as  the  ambient  group  and  to  the  set  S  as  the  target  set.  If  K  is  a  normal  subgroup 
of  G ,  then  Hf  =  G I  Kfp  is  a  group,  called  the  hidden  quotient  group,  and 
V  :G  — »  G  /  K  is  an  epimorphism,  called  the  hidden  epimorphism.  We  will  call  the 
above  diagram  the  hidden  subgroup  structure  of  the  map  tp:G  — >  S .  (See  [9,  59].) 


Remark.  The  underlying  intuition  motivating  this  formal  definition  is  as  follows:  Given 
a  natural  surjection  (or  epimorphism)  v  :G  ^Gl  K^,  a  "villain  with  malice  of 

forethought"  hides  the  algebraic  structure  of  v  by  intentionally  renaming  all  the  elements 
of  Gl K  ,  and  "tossing  in  for  good  measure"  some  extra  elements  to  form  a  set  S  and  a 

map  <p:G  — >  S  . 

The  hidden  subgroup  problem  can  be  stated  as  follows: 


1  By  saying  that  this  diagram  is  commutative,  we  m5ean  (p  =  l  °V  . 
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Hidden  Subgroup  Problem  (HSP).  Let  tp:G  — >  S  be  a  map  with  hidden  subgroup 
structure.  The  problem  of  determining  a  hidden  subgroup  of  G  is  called  a  hidden 

subgroup  problem  (HSP).  An  algorithm  solving  this  problem  is  called  a  hidden 
subgroup  algorithm. 

The  corresponding  quantum  form  of  this  HSP  is  stated  as  follows: 

Hidden  Subgroup  Problem  (Quantum  Version).  Let  cp:G  — >  S  be  a  map  with  hidden 
subgroup  structure.  Construct  a  quantum  implementation  of  the  map  (p  as  follows: 

Let  If  and  Hs  be  Hilbert  spaces  defined  respectively  by  the  orthonormal  bases 

j|  g) :  ge  G}  and  j|s) :  se  S}  and  let  s0  =  (p{f),  where  1  denotes  the  identity2  of  the 
ambient  group  A  .  Finally,  let  U  be  a  unitary  transformation  such  that 

hg®hs  >  hg®hs 

k)M  |^)|^>) 

Determine  the  hidden  subgroup  K  with  bounded  probability  of  error  by  making  as  few 
queries  as  possible  of  the  blackbox  U  v .  A  quantum  algorithm  solving  this  problem  is 
called  a  quantum  hidden  subgroup  (QHS)  algorithm. 


II.A.2.  The  generic  QHS  algorithm 

Let  (p :  G  — »  S  be  a  map  from  a  group  G  to  a  set  S  with  hidden  subgroup  structure. 
We  assume  that  all  representations  of  G  are  equivalent  to  unitary  representations3.  Let 

G  denote  a  complete  set  of  distinct  irreducible  unitary  representations  of  G .  Using 
additive  notation  for  G ,  we  let  1  denote  the  identity  of  G ,  and  let  s#  denote  its  image 

in  S  .  Finally,  let  1  denote  the  trivial  representation  of  G . 


Remark.  If  G  is  abelian,  then  G  becomes  the  dual  group  of  characters. 


The  generic  QHS  algorithm  is  given  below: 


2  We  are  using  multiplicative  notation  for  the  group  G  . 

3  This  is  true  for  all  finite  groups  as  well  as  for  a  large  class  of  infinite  groups. 
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Quantum  Subroutine  QRand  (#>) 


Step  0.  Initialization 


¥,)=  1)1 


Step  1.  Application  of  the  inverse  Fourier  transform  FG  *of  G  to  the  left  register 

vlGh 


geG 


where  lei  denotes  the  cardinality  of  the  group  G  . 


Step  2.  Application  of  the  unitary  transformation 

k. ) = -!=  E I W(g)> e » c  ® Hs 

yJ\G\  geG 


Step  3.  Application  of  the  Fourier  transform  FG  of  G  to  the  left  register 

1  (  ^1 

\¥,)  =  TFyi\7\Tmce  J>'(«)k}  \v(g))  =  ^Z\y\Trcice{\y)\<b(f)))^He®Hs, 


yeG  \g^G  J  |w  |  yeG 

where  |  y\  denotes  the  degree  of  the  representation  y ,  where  y'  denotes  the 

contragradient  representation  (i.e.,  y*  (g)  =  y^g'1 )  =  y(g)  ),  where 

M  M 


Trace ( / 1  y))  =  Z Z  U)|  Yv)>  and  where  ° ))  =  E  (tf)  |  P(tf))  • 


i=l  7=1 


geG 


Step  4.  Measurement  of  the  left  quantum  register  with  respect  to  the  orthonormal  basis 

{  I yu)-y^  G,l<i,j<\y\}. 

Thus,  with  probability 

\KJ\ 

the  resulting  measured  value  is  the  entry  yjj ,  and  the  quantum  system  "collapses" 
to  the  state 

M-ffr/)) 


k<)= 


S  Hd®Hs 


Step  5.  Output  y.j ,  and  stop. 
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II.A.3.  Shor’s  factoring  algorithm 

Briefly,  Shor’s  factoring  algorithm  is  described  as  follows: 

Let  N  be  the  integer  to  be  factored.  Choose  an  integer  Q  such  that  N2  <  Q  =  2*  <  IN2 , 
and  a  random  integer  t  which  with  high  probability  is  relatively  prime  to  N  .  Let 
<p :  Z  — >  Z N  be  the  map,  defined  by  the  n  t  mod  N  ,  from  the  additive  group  of 
integers  Z  into  the  integers  mod  N  (under  multiplication),  denoted  by  Z ;V  .  Shor’s 
algorithm  proceeds  to  find  the  hidden  subgroup  PZ  of  Z  of  all  multiples  of  P  ,  where 
P  is  the  unknown  period  of  the  map  (p  . 

This  is  accomplished  by  constructing  the  epimorphism  v :  Z  — »  ZQ  from  the  additive 
group  of  integers  Z  to  the  additive  group  of  integers  ZQ  modulo  Q  defined  by 
v:hb«  mod  Q ,  a  transversal4  i :  ZQ  — >  Z  for  the  epimorphism  v ,  and  a  map 
<p=  <p°i:ZQ  — >  Z  v which  is  a  “good  approximation”  to  the  map  (p . 

PZ  c  Z  — ^  ZN 


v  it  i  /*  (p  =  (p°i 

ZQ 

Then  Shor’s  algorithm  calls  the  quantum  subroutine  QRand($>)  to  produce  a  random 
rational  d!Q  which  is  with  high  probability  close  to  a  rational  of  the  form  y/P  .  By 
“d  IQ  close  to  y/P ,”  we  mean  that  d  /  Q  is  a  convergent  of  the  continued  fraction 
expansion  of  the  rational  y/P .  With  high  probability,  the  integers  y  and  P  are 
relatively  prime,  and  thus,  the  unknown  period  P  is  found. 

Remark.  The  mathematically  alert  reader  will  recognize  the  rational  d/Q  as  the 
character  xdtQ  :  ZQ  — >  R.  mod  1  of  ZQ  defined  by  n  mod  Q  I— >  e2mndlQ ,  and  the  rational 
y/P  as  the  character  %yjP  :  Zp  — >  R  mod  1  defined  by  n  modP  i-^  e2mnylp  _ 


4  By  a  transversal  of  V ,  we  mean  an  injection  l  such  that  l  °  V  =  idQ  is  the  identity  map  on  Z q  .  In 

other  words,  a  transversal  maps  each  element  of  Zg  into  an  element  of  the  coset  of  Z  representing  that 
coset. 
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II.A.4.  Wandering  Shor  algorithms,  a.k.a.,  vintage  shor  algorithms 

Wandering  Shor  algorithms  are  essentially  QHS  algorithms  on  free  abelian  finite  rank  n 
group  A  which,  with  each  iteration,  first  select  a  random  cyclic  direct  summand  Z  of 
the  group  A  ,  and  then  apply  one  iteration  of  the  standard  Shor  algorithm  to  produce  a 

random  character  of  the  “approximating”  finite  group  A-  Zg ,  called  a  group  probe5. 

Three  different  wandering  Shor  algorithms  are  created  in  [9],  The  first  two  wandering 
Shor  algorithms  given  in  [9]  are  quantum  algorithms  which  find  the  order  P  of  a 
maximal  cyclic  subgroup  of  the  hidden  quotient  group  H  .  The  third  computes  the  entire 

hidden  quotient  group  H . 

The  first  step  in  creating  a  wandering  Shor  algorithm  is  to  find  the  right  generalization  of 
the  transversal  i :  Zg  — >  Z  found  in  Shor’s  factoring  algorithm,  i.e.,  to  construct  the 

correct  transversal  from  Zg  to  a  free  abelian  group  A  of  rank  n .  For  this  reason, 
we  have  created  the  following  definition: 

Definition.  Let  A  be  the  free  abelian  group  of  rank  n,  let  Zg  be  the  cyclic  group  of 

order  Q ,  and  let  a  denote  a  chosen  generator  of  Z Q .  An  injection  i :  Zg  — >  A  is  said 
to  be  a  Shor  injection  provided  that: 

•  i^na'j  =  for  all  0  <  n  <  Q 

•  For  each  (free  abelian)  basis  ax  \a2 ,an '  of  A ,  the  coefficients 

A1,,A2,,...,Af  of  i{a)  =  Yjj^j'aj'  satisfy  gcd(V,V,...,/l„')  =  l. 

Remark.  Later,  when  we  construct  the  right  generalization  of  Shor  transversals  to  free 
groups  of  finite  rankn ,  we  will  see  that  a  Shor  transversal  is  nothing  more  than  a  2-sided 
Schreier  transversal.  Hence,  we  will  see  then  that  the  second  condition  of  the  above 

definition  simply  says  that  i  maps  the  generator  a  of  Zg  onto  a  generator  of  a  free 

direct  summand  Z  of  A  .  (For  more  details,  please  refer  to  section  II.A.9  of  this  report.) 
Remark.  In  [9],  we  show  how  to  use  the  extended  Euclidean  algorithm  to  construct  an 
epimorphism  v :  A  ZQ  having  i :  Zg  — >  A  as  its  transversal.  Consequently,  we  call  i 

a  Shor  transversal  for  the  epimorphism  v  . 


5  By  a  group  probe  A  ,  we  mean  an  epimorphic  image  of  the  ambient  group  A  . 
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Vlntage_Shor(  <p ,  Q,  n) 


Figure  1.  Flowchart  for  the  first  wandering  Shor  algorithm  (a.k.a.,  a  vintage  Shor 
algorithm).  This  algorithm  finds  the  order  P  of  a  maximal  cyclic  subgroup  of  the 

hidden  quotient  group  H . 

Flow  charts  for  the  three  wandering  Shor  algorithms  created  in  [9]  are  given  in  figures  1 
through  3.  In  [9],  these  were  also  called  vintage  Shor  algorithms. 
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A.ltl_  Vintage _ Shor(  <p  ,  Q,  n,  K) 


Figure  2.  Flowchart  for  the  second  wandering  Shor  algorithm  (a.k.a.,  a  vintage 
Shor  algorithm).  This  algorithm  finds  the  order  P  of  a  maximal  cyclic  subgroup  of 

the  hidden  quotient  group  H  . 


The  algorithmic  complexities  of  the  above  wandering  Shor  algorithms  is  given  in  [9], 
For  example,  the  first  wandering  Shor  algorithm  is  of  time  complexity 

0(n2(lg7V)3(lglg7V)"+1), 

where  n  is  the  rank  of  the  free  abelian  group  A  .  This  can  be  readily  deduced  from  the 
following  abbreviated  flowchart  given  in  figure  4. 
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Alt2_Vintage_Shor(  <p ,  Q,  n  ,  K) 


Figure  3.  Flowchart  for  the  third  wandering  Shor  algorithm,  a.k.a.,  a  vintage  Shor 
algorithm.  This  algorithm  finds  the  entire  hidden  quotient  group  H  . 


jstep  lj  O(n) 


Figure  4.  Abbreviated  flowchart  for  the  first  wandering  Shor  algorithm. 
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II.A.5.  A  continuous  (variable)  Shor  algorithm 


A  continuous  variable  Shor  algorithm  was  created  in  [13,  24],  By  a  continuous 
variable  Shor  algorithm,  we  mean  a  quantum  hidden  subgroup  algorithm  that  finds  the 
hidden  period  P  of  an  admissible  function  (p :  R  — >  R  from  the  reals  R  to  itself. 

Remark.  By  an  admissible  function,  we  mean  a  function  belonging  to  any  sufficiently 
well  behaved  class  of  functions.  For  example,  the  class  of  functions  which  are  Lebesgue 
integrable  on  every  closed  interval  of  R .  There  are  many  other  classes  of  functions  that 
work  equally  as  well. 

Actually,  the  papers  [13,  24]  give  in  succession  three  such  continuous  Shor  algorithms, 
each  successively  more  general  than  the  previous. 

For  the  first  algorithm,  we  assume  that  the  hidden  period  P  is  an  integer.  The  algorithm 
is  then  constructed  by  using  rigged  Hilbert  spaces,  linear  combinations  of  Dirac  delta 
functions,  and  a  subtle  extension  of  the  Fourier  transform  in  the  generic  QHS  subroutine 
QRand(#>)  ,  which  has  been  described  previously  in  section  11.A.2  of  this  report.  In  Step 
5  of  this  algorithm,  the  observable 


is  measured,  where  Q  is  an  integer  chosen  so  that  Q  >  IP1 .  ft  then  follows  that  the 


output  of  this  algorithm  is  a  rational  ml  Q  which  is  a  convergent  of  the  continued 
fraction  expansion  of  a  rational  of  the  form  n  /  P . 

The  above  quantum  algorithm  is  then  extended  to  a  quantum  algorithm  (the  second 
algorithm)  that  finds  the  hidden  period  P  of  functions  (p :  R  — >  R ,  when  Pisa  rational. 

Finally,  the  second  algorithm  is  extended  to  the  third  algorithm  which  finds  the  hidden 
period  P  of  functions  (p :  R  — »  R ,  whenP  is  an  arbitrary  real  number.  We  point  out  for 
the  third  and  last  algorithm  to  work,  we  must  impose  a  very  restrictive  condition  on  the 
map  tp :  R  — >  R ,  i.e.,  the  condition  that  the  map  (p  is  continuous. 

II.A.6.  Three  techniques  for  creating  new  QHS  algorithms,  i.e.,  lifting,  pushing, 
and  duality 

We  have  already  seen  one  ingredient  of  Shor’s  factoring  algorithm  that  we  were  able  to 
generalize  to  create  new  QHS  algorithms.  Specifically,  we  have  generalized  the 
transversal  i:  ZQ  — >  Z  found  in  Shor’s  original  algorithm  to  a  transversal  i :  Zg  — >  A 

(called  a  Shor  transversal)  from  the  finite  cyclic  group  Zeto  the  free  abelian  group  A 
of  finite  rank.  This  enabled  us  to  create  new  quantum  algorithms,  called  wandering  Shor 
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algorithms.  (See  section  II.A.4.)  Later,  in  section  II.A.9,  we  will  see  how  this  ingredient 
can  be  further  generalized  to  a  2-sided  Schreier  transversal,  thereby  enabling  us  to 
create  new  QHS  algorithms  on  free  groups. 

In  this  section,  we  seek  to  generalize  other  ingredients  found  within  Shor’s  original 
algorithm  with  the  intention  of  using  these  generalizations  to  create  even  more  new 
quantum  algorithms.  For  that  reason,  we  give  the  following  definitions: 

Definition.  A  map  tp\G  — »  S  from  a  group  G  to  a  set  S  is  said  to  be  the  push  of  a  map 
<p:G  — >  S  from  a  group  G  to  S ,  written 

(p  -  Push[(p) 

if  there  exists  an  epimorphism  v  :G  —>  G  of  G  onto  G  and  a  transversal  i:G  — >  (7  of  v 
such  that  <p  =  <p°i . 


Thus,  in  terms  of  the  above  definition,  we  can  now  see  that  Shor’s  algorithm,  as 
described  in  section  H.A.3.,  is  an  algorithm  created  by  pushing  the  hidden  subgroup 

problem  (HSP)  tp :  Z  — »  S  to  Push(<p )  =  <p:ZQ  ^  ZN  ,  and  then  calling  QRand  (  j  to 
produce  a  character  %  IQ  of  the  group  Zg .  The  character  %yjQ  is  with  high  probability 
close  to  a  primitive  character  %djP  of  the  hidden  subgroup  Zp .  Hence,  the  character 
Xd,P  can  be  found  with  the  standard  continued  fraction  algorithm. 


We  next  note  that  the  above  generic  definition  of  pushing  suggests  a  second  procedure 
which  can  in  turn  also  be  used  for  creating  even  more  QHS  algorithms,  namely: 

Definition.  A  map  cp\G  — »  S  from  a  group  G  to  a  set  S  is  said  to  be  the  lift  of  a  map 
(p\G  — >  S  from  a  group  G  to  S ,  written 

<P  =  Lift(<p) 

if  there  exists  a  morphism  jU:G  — >  (7  of  G  into  G  such  that  (p  =  (p°  p . 


There  is  yet  a  third  generic  procedure  that  we  have  used  to  create  new  QHS  algorithms  — 
duality. 

Let  G  be  an  abelian  group,  and  let  G  denote  its  dual  group  of  characters.  Then  if  we 
are  creating  QHS  algorithms  for  hidden  subgroup  problems  (HSPs)  of  the  form 
tp'.G  ^  S ,  we  might  as  well  use  the  same  procedures  to  create  QHS  algorithms  for 

HSPs  of  the  form  (p\G  — >  S  .  We  refer  to  this  method  of  creating  new  QHS  algorithms 

as  duality. 
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In  the  sections  to  follow,  we  will  show  how  lifting,  pushing,  and  duality  are  used  to 
create  new  algorithms.  In  particular,  in  section  II.A.7,  we  use  lifting  and  duality  to  create 
three  new  QHS  algorithms,  namely,  the  lifted  Shor  algorithm,  the  circle  algorithm,  and 
the  dual  Shor  algorithm.  In  section  II.A.9,  we  will  use  pushing  and  lifting  to  create  QHS 
algorithms  on  free  groups. 


II.A.7.  The  lifted  Shor,  the  quantum  circle,  and  the  dual  Shor  algorithms 

The  methods  of  lifting  and  duality  given  in  the  previous  section  are  used  in  [13,  19,  24]  to 
create  the  following  QHS  algorithms: 

•  The  lifted  Shor  algorithm  —  Created  from  the  Shor  algorithm  by  lifting 

•  The  circle  algorithm  —  Created  from  the  lifted  Shor  algorithm  by  duality 

•  The  dual  Shor  algorithm  —  Created  from  the  circle  algorithm  by  lifting 


A  roadmap  for  creating  these  algorithms  is  given  in  figure  5. 


Figure  5.  Roadmap  for  creating  new  QHS  algorithms. 

As  its  name  suggests,  the  lifted  Shor  algorithm  is  created  by  lifting  the  hidden  subgroup 
problem  (HSP)  :  Zg  — >  Z^  to  a  HSP  (p :  Z  — >  S ,  thereby  producing  a  QHS  algorithm 

which  is  essentially  a  “distillation”  of  Shor’s  original  algorithm.  Next,  duality  is  used  to 
create  the  quantum  circle  algorithm  by  devising  a  QHS  algorithm  for  the  HSP 
(p :  M  /  Z  — >  S  on  the  dual  group  R/Z  of  the  additive  group  of  integers  Z .  (By  R  /  Z , 
we  mean  the  additive  group  of  reals  mod  1,  which  is  isomorphic  to  the  multiplicative 
group  | e,e  :  0  <  6  <  2#} ,  i.e.,  the  unit  circle  in  the  complex  plane.)  Finally,  the  dual 

Shor  algorithm  is  created  by  lifting  the  HSP  (p :  R/Z  — >  S  to  the  HSP  <p:ZQ^S . 
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It  is  not  clear  whether  or  not  the  lifted  Shor  algorithm  and  the  quantum  circle  algorithm6 
are  physically  implementable  quantum  algorithms.  Even  so,  they  play  an  important  role 
as  two  essential  stepping  stones  to  the  very  implementable  dual  Shor  algorithm. 

For  detailed  descriptions  of  each  of  these  quantum  algorithms,  i.e.,  the  lifted  Shor,  the 
quantum  circle,  and  the  dual  Shor  algorithms,  the  reader  is  referred  to  [13,  19,  24], 

We  give  below  brief  descriptions  of  the  quantum  circle  and  the  dual  Shor  algorithms. 

For  the  quantum  circle  algorithm,  we  make  use  of  the  following  spaces  (each  of  which 
is  used  in  quantum  optics) 

•  The  rigged  Hilbert  space  HRIZ  with  orthonormal  basis  j|  x) :  x  e  R.  /  Z} .  By 

“orthonormal”  we  mean  that  (jc  |  =  <£(  x  -  j) ,  where  “  8  ”  denotes  the  Dirac 

delta  function.  The  elements  of  HR/Z  are  formal  integrals  of  the  form 

/(jc)|  x) .  (The  physicist  Dirac  in  his  classic  book  on  quantum  mechanics 

refers  to  these  integrals  as  infinite  sums.) 

•  The  complex  vector  space  Hz  of  formal  sums 

c  V«e  zj 

with  orthonormal  basis  {|  w) :  n  e  Z} .  By  “orthonormal”  we  mean  that 
(n\m)  =  8nm  ,  where  Snm  denotes  the  Kronecker  delta. 


We  will  now  design  an  algorithm  which  solves  the  following  hidden  subgroup  problem: 

Hidden  Subgroup  Problem  for  the  Circle.  Let  (p :  R  /  Z  — >  C  be  an  admissible 
function  from  the  circle  group  M/Z  to  the  complex  numbers  C  with  hidden  rational 
period  ere  Q/Z,  where  a^QI'L  denotes  the  rational  circle,  i.e.,  the  rationals  mod  1. 

Remark.  By  an  admissible  function,  we  mean  a  function  belonging  to  any  sufficiently 
well  behaved  class  of  functions.  For  example,  the  class  of  functions  which  are  Lebesgue 
integrable  on  M  /  Z .  There  are  many  other  classes  of  functions  that  work  equally  as  well. 

Proposition.  Let  a- ax!  a2  (with  gcd(«j,a2)  =  1)  be  a  period  of  a  function 

(p :  M/Z  — >  C .  Then  1/  a2  is  also  a  period  of  <p .  Hence,  the  minimal  rational  period  of 

(p  is  always  a  reciprocal  integer  mod  1 . 


6  There  is  a  possibility  that  the  quantum  circle  algorithm  may  have  a  physical  implementation  in  terms  of 
quantum  optics. 
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The  following  quantum  algorithm  finds  the  reciprocal  integer  period  of  the  function  (p . 


Circle- Algorithm(  (p ) 


Step  0.  Initialization 


Step  1.  Application  of  the  inverse  Fourier  transform  F  1  ®  1 

| y, )  =  | dx  e2’1- 1  x) 1 0>  =  |  dx  \  x) | 0)  <=  ®  Hc 

Step  2.  Application  of  the  unitary  transformation  U9  :  |  x)  |  u)  I— >  \  x)  \  u  +  (p{x)) 

I  Yi ) —  ^  I  *)  I  ^A))  ^  -®r/z  ®  Hc 
Step  3.  Application  of  the  Fourier  transform  F  <8>  1 


k*>=Z^  e~2ninx  I n)  I  ?**)) = E I  ^  e~21tinx  I  ^x)) e  © 


H 


neZ 


neZ 


m 

Remark.  Letting  xm  —  x - ,  we  have 

a 

a- 1  (m+l)/a 

|  dx  eir,nx  |  <p( x))  =  X  {  dx  e-lKinx  |  ^(jc)) 

m= 0  ml  a 


Tin[Xm  +  tj] 


1La  -2nin\  xm  +—  | 

£ 

w=0  o 

/  «-i  A  i/« 


x„.  +- 


in 


Yje~lKinm,a  J  dx e~2Kinx\(p(x)) 


\m=0 


where  1  /  a  is  the  unknown  reciprocal  period.  But 


£ 


a- 1 

< 

m=0 


-Ininm  /  a 


= 


//=()  mod  a 


\  a  if  n-  0  mod  a 
0  otherwise 


Hence, 


\¥3)  =  Yd\n)^dxe~1,l,HX\<p{x))  =  ZlwK=0mod«  I  j  dx  e-2*inx\(p(x)) 

veZ  \ntZ  J  o 

'  \flla  A 

^|fa)  J  dx  e~2mnx\(p(x)^  -  ^|la)|fl(lfl)) 
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Step  4.  Measurement  of 


with  respect  to  the  observable 


Ysn\n)(n 


to  produce  a  random  eigenvalue  £a  . 

Remark.  The  above  quantum  circle  algorithm  can  be  extended  to  a  quantum  algorithm 
which  finds  the  hidden  period  P  of  a  function  (p :  R  /  Z  — »  C  ,  when  P  is  an  arbitrary 
real  number  mod  1 .  But  in  creating  this  extended  quantum  algorithm,  we  must  impose  a 
very  restrictive  condition  on  the  map  (p :  R  /  Z  — »  C ,  i.e.,  the  condition  that  the  map  (p  is 
continuous. 

We  now  give  a  brief  description  of  the  dual  Shor  algorithm. 

The  dual  Shor  algorithm  is  a  QHS  algorithm  created  by  making  a  discrete  approximation 
of  the  quantum  circle  algorithm.  More  specifically,  it  is  created  by  lifting  the  QHS  circle 
algorithm  for  cp :  R  /  Z  — »  C  to  the  finite  cyclic  group  Zg ,  as  illustrated  in  the 

commutative  diagram  given  below: 


jui  \  (p-  Push((p)  =  (po  ju 


R/Z  — 5 


Intuitively,  we  “approximate”  the  circle  group  R/Z  with  the  finite  cyclic  group  Zg  by 
identifying  Zg  with  the  following  additive  group 


and  by  identifying  the  hidden  subgroup  Zp  with  the  additive  group 


where  P  =  a2. 

This  is  a  physically  implementable  quantum  algorithm.  It  is  actually  faster  than  Shor’s 
algorithm.  For  the  last  step  of  Shor’s  algorithm  uses  the  standard  continued  fraction 
algorithm  to  determine  the  unknown  period.  On  the  other  hand,  the  last  step  of  the  dual 
Shor  algorithm  uses  the  much  faster  Euclidean  algorithm  to  compute  the  greatest 
common  divisor  of  the  integers  £x  a,  i2a,  l3a, ...  ,  thereby  determining  the  desired 

reciprocal  integer  period  1/ a .  For  more  details,  please  refer  to  [13,  19,  14]. 
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II.A.8.  A  QHS  algorithm  for  Feynman  integrals 

We  now  discuss  a  QHS  algorithm  based  on  Feynman  path  integrals.  This  quantum 
algorithm  was  developed  at  the  Mathematical  Sciences  Research  Institute  (MSRI)  in 
Berkeley,  California  when  the  PI  was  challenged  by  being  invited  to  give  a  talk  on 
Feynmann  path  integrals  and  quantum  computing  at  an  MSRI  conference  on  Feynman 
path  integrals. 

Until  recently,  the  PI  thought  that  the  quantum  algorithm  to  be  described  below  was  a 
highly  speculative  quantum  algorithm,  because  the  existence  of  Feynman  path  integrals  is 
very  difficult  (if  not  impossible)  to  determine  in  a  mathematically  rigorous  fashion.  But 
much  to  the  Pi’s  surprise,  Jeremy  Becnel  in  his  doctoral  dissertation  [36]  has  succeeded 
in  creating  a  firm  mathematical  foundation  for  this  algorithm. 

We  should  mention,  however,  that  the  physical  implementability  of  this  algorithm  is  still 
yet  to  be  determined. 

Definition.  Let  Paths  be  the  real  vector  space  of  all  continuous  paths  x  = [0,1]  ->  R" 
which  are  Is  with  respect  to  the  inner  product 

x»y  =  |  ds  x(s)y(s) 

with  scalar  multiplication  and  vector  sum  defined  as 

•  (Ax)(s)  =  Ax(s) 

•  (x  +  y)(s)  =  x(s)  +  y(s) 


We  wish  to  create  a  QHS  algorithm  for  the  following  hidden  subgroup  problem: 

Hidden  Subgroup  Problem  for  Paths .  Let  (p :  Paths  — »  C  be  a  functional  with  a 
hidden  subspace  V  of  Paths  such  that 

tp(x  +  v)  =  <p(x)  VveF 


Our  objective  is  to  create  a  QHS  algorithm  which  solves  the  above  problem,  i.e.,  which 
finds  the  hidden  subspace  V . 

Definition.  Let  HPaths  be  the  rigged  Hilbert  space  with  orthonormal  basis 
j|x):xe  Paths},  and  with  bracket  product  (x|y)  =  <?(x-y). 

We  will  use  the  following  observation  to  create  the  QHS  algorithm: 


18 


Observation.  Paths  -  [J  K(v  +  ^/±)  >  where  V1  denotes  the  orthogonal  complement  of 
the  hidden  vector  subspace  V . 


The  QHS  algorithm  for  Feynman  path  integral  is  now  given  below: 

Feynman  ($>) 


Step  0.  Initialize 


Step  1.  Apply  F  ®  1 


|r,)=  J  Dx  e"M|x)|0)  =  J  Dx  |x)|0 


Paths 


Paths 


Step  2.  Apply  U  :  |jc)|m)  i— >  |jc)|«  +  tp(x)) 

| Wi)=  J  -®*|*)|p(*)) 


Paths 


Step  3.  Apply  F  ®  1 


tyi)  =  \  Dy  \  Dxe-**->\y)\<p(x) 

Paths  Paths 

=  {  Dy\y)  {  Dxe^^x) 


Paths  Paths 


But 


However, 


{  Dxe-1*ix-,,\(p(x))  =  \Dv  J  Dxe-2™-y\(p(x ) 

Paths  V  v+VL 

=  \Dv  |  Dx e-2m(v+xhy  |^?(v  +  jc) 

V  V1 

=  \Dv  e-lKi"‘y  J  Dx  e~2,cix'y\(p{x) 

V  V1 

\Dve~2™'y  =\  DuS(y-u) 


So, 
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k,)=  J  Dy\y)\Dve-,’‘-"\DXe-M-’\<p(x)) 

Pathsn  V  V1 

=  {  Dy\y)  ^Dud(y-u)  ^  Dx  e~2nix‘y\(p(x)) 

Pathsn  V1  V1 

=  |  Z)ujuj  j  Dx e~2mx'u\(p[x)^j 

vL  vL 

-  |  Du  |*f)|£2(n)) 

v1 

Step  4.  Measure 

\v3)=  J 

V1 

with  respect  to  the  observable 

A  =  |  Dn  w|w)(w| 

Paths 

to  produce  a  random  element  of  V1 


The  above  algorithm  suggests  an  intriguing  question.  Can  the  above  QHS  Feynman 
integral  algorithm  be  modified  in  such  a  way  as  to  create  a  quantum  algorithm  for  the 
Jones  polynomial?  In  other  words,  can  it  be  modified  by  replacing  Paths  with  the  space 
of  gauge  connections,  and  making  suitable  modifications? 

This  question  is  motivated  by  the  fact  that  the  integral  over  gauge  transformation 

y/{K)  =  \DAy/{A)wK{A) 
looks  very  much  like  a  Fourier  transform,  where 

WK(A)  =  tr[p™v{§KA)) 

denotes  the  Wilson  loop  over  the  knot  K  . 


II.A.9.  Non-abelian  QHS  algorithms  --  A  simplification 

Let  (p :  G  — »  S  be  a  map  with  hidden  subgroup  structure  from  a  finitely  generated  (f.g.) 
non-abelian  group  G  to  a  set  S .  We  assume  that  the  hidden  subgroup  A  is  a  normal 
subgroup  of  G . 
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Under  the  above  assumptions,  we  now  outline  why  one  need  only  study  non-abelian  QHS 
algorithms  on  free  groups.  More  specifically,  we  briefly  describe  the  following  results: 

•  Every  hidden  subgroup  problem  tp :  G  — >  S  on  a  non-abelian  f.g.  group  G  can  be 

lifted  to  a  hidden  subgroup  problem  (p:  F  — >  S  on  a  f.g.  free  group  F ,  and  in 
turn 

•  Every  polytime  QHS  algorithm  for  (p:F  — >  S  can  be  transformed  via  the 
Reidemeister-Schreier  theorem  into  a  polytime  QHS  algorithm  for  the  original 
hidden  subgroup  problem  (p:G  — >  S  . 

But  what  do  we  mean  by  a  free  group? 

(Universal)  Definition.  A  f.g.  group  F  is  said  to  be  free  if  there  exists  a  set  of 
generators  X  =  {x1,x2,...,xn}  such  that,  for  every  group  G  and  for  every  map 
f  :X  — >  (7  of  the  set  X  into  the  group  G ,  the  map  f  extends  to  a  morphism 
f  :  F  — >  G  .  We  call  the  set  X  a  free  basis  of  the  group  F ,  and  frequently  denote  the 
group  F  by  F  (x1,x2,...,xn) .  It  follows  from  this  definition  that  the  morphism  f  is 
unique. 


The  intuitive  idea  encapsulated  by  this  definition  is  that  a  free  group  is  an  unconstrained 
group  (very  much  analogous  to  a  physical  system  without  boundary  conditions.)  In 
other  words,  a  group  is  free  provided  it  has  a  set  of  generators  such  that  the  only  relations 
among  those  generators  are  those  required  for  F  to  be  a  group.  For  example, 

•  x;x f1  =1  is  an  allowed  relation 

•  x.x,  -  x.x.  is  not  an  allowed  relation  for  i  ^  / 

1  J  J  1  J 

•  xf3  =  1  is  not  an  allowed  relation 


Definition.  Let  F  ( x2 ,  x2 , . . . ,  xn )  be  a  free  group  with  free  basis  xl,x2,...,xn.  Then  a 
word  is  a  finite  string  of  the  symbols  x,,  x,_l,  x2,  xf1, ...  ,xn,  xfl .  A  reduced  word  is  a 
word  in  which  there  is  no  substring  of  the  form  XjXf1  or  xf'xj .  Two  words  are  said  to 
be  equivalent  if  one  can  be  transformed  into  the  other  by  applying  a  finite  number  of 
substring  insertions  or  deletions  of  the  form  XjXf1  or  xf'xj .  We  denote  an  arbitrary 

word  w  by  w  -  ala2---ae  ,  where  each  a.  -  xk±l  . 

For  example,  x2x1“1x1x1“1x5“1x5“1x5“1x5  is  a  word  which  is  equivalent  to  the  reduced 
word  x2x1“1x5"1x5'i . 
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It  easily  follows  that  a  free  group  F  ( x1 ,  x2 , . . . ,  xn )  is  nothing  more  than  the  set  of 

reduced  words  together  with  the  obvious  definition  of  product,  i.e.,  concatenation  with 
full  reduction. 

Definition.  Let  G  be  a  group.  A  group  presentation 

(x1,x2,...,xn:rl,r2,...,rm) 

for  G  is  a  set  of  free  generators  x1,x2,...,xn  of  a  free  group  F  and  a  set  of  words 
rx,r2,...,rn  in  F  ( xx ,  x2 , . . . ,  xn ) ,  called  relators,  such  that  the  group  G  is  isomorphic  to 
the  quotient  group  F(x1,x2,...,xn)/Cons(r1,r2,...,rit),  where  Cons  ( /' ,  r2 , . . . ,  rn )  , 
called  the  consequence  of  rvr2,...,rn>  is  the  smallest  normal  subgroup  of 
F(xx, x2,...,xn)  containing  the  relators  rx,r2,...,rn. 


The  intuition  captured  by  the  above  definition  is  that  x1,x2,...,xn  are  the  generators  of 
G ,  and  rx  =  1,  r2  =  1, . . . ,  rn  =  1  is  a  complete  set  of  relations  among  these  generators,  i.e., 
every  relation  among  the  generators  of  G  is  a  consequence  (derivable  from)  the  relations 
rx  =  1,  r2  =  1, . . . ,  rn  =  1 .  For  example, 

•  (xx,x2,...,xn :)  and  (xx,x2,...,xn  :  x^-^xfxf^x^xf'xf')  are  both 
presentations  of  the  free  group  F  ( xx ,  x2 , . . . ,  xn ) 

•  ( x :  xQ )  and  ( x :  x“ ,  xh  j  are  both  presentations  of  the  cyclic  group  ZQ  of  order 
Q ,  where  a  and  b  are  integers  such  that  gcd  («,/>)  =  Q . 

•  |xj,x2 :  Xj3,x22,(xjX2)2  j  is  a  presentation  of  the  symmetric  group  S3  on  three 
symbols. 

The  procedure  we  are  about  to  describe  makes  use  of  the  Reidemeister-Schreier  (R-M) 
theorem,  a  well  known  theorem  in  combinatorial  group  theory.  Given  a  presentation 
(xx,x2,...,xn  :rx,r2,...,rm)  of  a  group  G  andasubgroup  K  of  G ,  the  R-M  theorem 

gives  an  algorithm  for  computing  a  presentation  of  the  subgroup  K  .  This  R-M  theorem 
algorithm  is  constructed  using  a  Schreier  transversal. 

But  what  is  a  Schreier  transversal? 

Definition.  A  set  W  of  reduced  words  in  a  free  group  F  =  F  ( xx ,  x2 , . . . ,  xn )  is  said  to 

be  a  2-sided  Schreier  system  provided 

•  w  =  axa2  ■  ■  •al_lae  e  W  =>  wLeft  =  axa2  ■  -a^  e  W ,  and 

•  w  =  axa2  •  •  •  ae_xat  e  W  =>  wRight  =a2  -  a^a,  e  W 
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Given  an  epimorphism  v  :  F  — >  (7  of  the  free  group  F  onto  a  group  G ,  a  2-sided 
Schreier  transversal  T  :G  — »  F  for  v  is  a  transversal  for  v  for  which  there  exists  a  2- 
sided  Schreier  system  such  that  t(G)  =  W . 

Remark.  Most  surprisingly,  the  concept  of  a  2-sided  Schreier  transversal  provides  a 
very  natural  way  of  extending  Shor  ’s  original  factoring  algorithm  to  a 
quantum  algorithm  on  free  groups.  For  the  transversal  found  in  Shor 's  original  quantum 
algorithm  and  in  the  wandering  Shor's  algorithm  are  nothing  more  than  abelian  2-sided 
Schreier  transversals!  (See  sections  II.A.3  and II.A.4.) 


Finally,  we  are  now  prepared  to  describe  the  central  idea  of  this  section. 

Given  the  hidden  subgroup  problem  (HSP)  <p:G  — »  S  with  hidden  subgroup  K ,  we 

construct  a  QHS  algorithm  for  (p  using  the  following  procedure: 

STEP  1 .  Construct  a  presentation  ( x, , . . . ,  xn :  r, , . . . ,  rm )  for  the  group  G  ,  and  let 
v :  F  ( Xj ,  x2 , . . . ,  xm )  — »  G  denote  the  epimorphism  corresponding  to  this 
presentation. 

STEP  2.  Construct  the  lifted  HSP  (p  =  Lift {(p)  =  (p°v :  F  — »  S  .  (See  section  II. A.6.) 

STEP  3.  Use  QRand^cp^  to  find  the  hidden  subgroup  K  =  v_1  ( A)  of  the  HSP  (p . 

(See  section  II.A.7.) 

STEP  4.  Construct  a  2-sided  Schreier  transversal  t:G  — »  F  for  v . 

STEP  5.  With  the  above  transversal  r ,  use  the  R-M  Theorem  to  construct  a  presentation 
of  the  hidden  subgroup  K  of  the  original  HSP  (p . 


More  details  about  this  approach  to  creating  new  non-abelian  QHS  algorithms  can  be 
found  in  the  forthcoming  paper  [32], 


II.A.10.  Is  Grover’s  algorithm  a  QHS  algorithm? 

Is  Grover’s  algorithm  a  quantum  hidden  subgroup  algorithm?  Does  Grover’s  algorithm 
have  some  symmetries  that  we  can  exploit? 
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The  problem  solved  by  Grover’s  algorithm  [7,  51,  52,  53]  is  that  of  finding  an  unknown 
integer  label  j0  in  an  unstructured  database  with  items  labeled  by  the  integers: 

0,1,2,  — 1  =  2"— 1, 


given  the  oracle 


/(/)  = 


if  j  =  j0 

otherwise 


Let  H  be  the  Hilbert  space  with  orthonormal  basis  1 0) ,  1 1) ,  1 2) , . . . ,  |  N  - 1) .  Grover ’  s 
oracle  is  essentially  given  by  the  unitary  transformation 

/  :  H  H 

\h) 

I i)  «  (-lfW 

where  /|4)=/-2Ia)UI  is  inversion  in  the  hyperplane  orthogonal  to  |  j)  .  Let 

W  denote  the  Hadamard  transformation  on  the  Hilbert  space// .  Then  Grover’s 
algorithm  is  as  follows: 

STEP  0.  (Initialization) 

k  ^0 

STEP  1.  Loop  until  k  ~  7t4n  /  4 

k)  <-  ek)=-W|„pU)k> 

k  ^ —  k+1 

STEP  2.  Measure  |  y/)  with  respect  to  the  standard  basis 
|0),|1),|2),...,|7V-1) 
to  obtain  the  unknown  state  y(l)  with 

Prob  > 1 -  — 

N 

But  where  is  the  hidden  symmetry  in  Grover’s  algorithm? 

Let  SN  be  the  symmetric  group  on  the  symbols  0, 1, 2, . . .,  TV  - 1 .  Then  Grover’s 
algorithm  is  invariant  under  the  hidden  subgroup  Stabh  -  jg €  SN  :  g(j0)  =  jt) }  SN , 
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called  the  stabilizer  subgroup  for  j0 ,  i.e.,  Grover’s  algorithm  is  invariant  under  the 
group  action 

Stub-  xH  — ^  H 

J  0 

Moreover,  if  we  know  the  hidden  subgroup  Stabj  ,  then  we  know  j0 ,  and  vice  versa.  In 
other  words,  the  problem  of  finding  the  unknown  label  j0  is  informationally  the  same  as 
the  problem  of  finding  the  hidden  subgroup  Stabh . 

Let  (iy)e  SN  denote  the  permutation  that  interchanges  integers  i  and  j ,  and  leaves  all 
other  integers  fixed.  Thus,  (//’)  is  a  transposition  if  /  ^  j ,  and  the  identity  permutation  1 
if  i  =  j. 

Proposition.  The  set  {  (Oy0),(ly0),(2y#),...,((/V-l)y#)  }  is  a  complete  set  of 
distinct  coset  representatives  for  the  hidden  subgroup  Stab •.  of  SN,  i.e.,  the  coset  space 
SN  /  Stabj  is  given  by  the  following  complete  set  of  distinct  cosets: 

SNl Stab.  =  {  ( 0  f ) Stab. , ( 1  f ) Stab. , ( 2  f ) Stab. , . . . , ( ( /V - 1 ) y# ) Stab.  } 


We  can  now  see  that  Grover’s  algorithm  is  a  hidden  subgroup  algorithm  in  the  sense  that 
it  is  a  quantum  algorithm  which  solves  the  following  hidden  subgroup  problem: 

Grover’s  Hidden  Subgroup  Problem.  Let  tp:  SN  S  be  a  map  from  the  symmetric 
group  SN  to  a  set  S  =  {  0,1,2,...,  iV  — 1}  with  hidden  subgroup  structure  given  by  the 
commutative  diagram 

SN  -*->  S 

v,  \  , 

Jo 

S,,  /  Stab. 

Jo 

where  v.  :  SN  — >  SN  /  Stab  ,  is  the  natural  surjection  of  SN  on  to  the  coset  space 
S x.  /  Stab.  ,  and  where 

Jo 

i  •  S  N  /  Stub .  — )  S 

iy  Jo 

(j  jo) Stab h  i->  j 

is  the  unknown  relabeling  (bijection)  of  the  coset  space  SN  /  Stab,  onto  the  set  S .  Find 
the  hidden  subgroup  Stab .  with  bounded  probability  of  error. 
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Now  let  us  compare  Shor's  algorithm  with  Grover's. 

From  section  II.A.3,  we  know  that  Shor’s  algorithm  [5,  9,  71,  72]  solves  the  hidden 
subgroup  problem  (p :  Z  — >  ZN  with  hidden  subgroup  structure 

Z  — ^  ZN 

v\  Si 

7LIP7L 

Moreover,  as  stated  in  section  II.A.3,  Shor  creates  his  algorithm  by  pushing7  the  above 
hidden  subgroup  problem  :  Z  — >  Z^  to  the  hidden  subgroup  problem 

(p :  Zg  — >  Z  v  (called  Shor’s  oracle),  where  the  hidden  subgroup  structure  of  (p  is  given 
by  the  commutative  diagram 

Z  — ^  ZN 

a  \\  t  S(p  =  (p°r  , 

Z  Q 

where  a  is  the  natural  epimorphism  of  Z  onto  Zg ,  and  where  t  is  Shor’s  chosen 
transversal  for  the  epimorphism  a . 


Surprisingly,  Grover's  algorithm,  viewed  as  an  algorithm  that  solves  the  Grover  hidden 
subgroup  problem,  is  very  similar  to  Shor’s  algorithm. 


Like  Shor’s  algorithm,  Grover’s  algorithm  solves  a  hidden  subgroup  problem,  i.e.,  the 
Grover  hidden  subgroup  problem  (p\SN^S  with  hidden  subgroup  structure 


V\  Si 


SN I  Stab. 

iy  J  0 

where  S  =  {0,1,2,...,  A-l}  denotes  the  set  resulting  from  an  unknown  relabeling 
(bijection) 

(j  j*)  Stab hv*j 


of  the  coset  space 

SNIStabh  ={  ( 0  y(1 )  Stab. ,  ( 1  ,/# )  Stab. ,  ( 2  y(l )  Stab. , . . . ,  ( ( TV  - 1  ),/# )  Stab.  }. 


Also,  like  Shor’s  algorithm,  we  can  think  of  Grover’s  algorithm  as  one  created  by 
pushing  the  Grover  hidden  subgroup  problem  (p :  SN  — »  S  to  the  hidden  subgroup 


7  See  Section  II.A.6  for  a  definition  of  pushing. 


26 


problem  (p :  SN  /  Stab  .  — »  5 ,  where  the  pushing  is  defined  by  the  following  commutative 


diagram 

— 

S  -  SN  / Stab. 

iy  J  0 

a  \\  t 

o 

II 

\  §- 
\ 

SN  /  Stab, 

where  a:SN  — >  5^  / denotes  the  natural  surjection  of  5^ onto  the  coset  space 
Av  /  .SVaA,, ,  and  where  t:Sn/  Stab,  — >  denotes  the  transversal  of  or  given  by 

iSyy  /  Stab,  — > 

(7  0)Stab,  (y  0)' 


Again  also  like  Shor’s  algorithm,  the  map  (p  given  by 

Sy  /  Stab,  — >  Sy  /  Stab  j  —  iS1 

(y  0) Stab,  (y  y0 )  Stab.  =  j 

is  (if  y#  ^  0 )  actually  a  disguised  Grover’s  oracle.  For  the  map  (p  can  easily  be  shown  to 
simply  to 


tp((j  0)Stab0)  = 


(y  (i)Stab. 
Stab. 


ifj  =  J. 

otherwise 


which  is  informationally  the  same  as  Grover’s  oracle 

f  j  if  J  =  ./'o 
1 1  otherwise 


/(/)= 


Hence,  we  can  conclude  that  Grover's  algorithm  is  a  quantum  algorithm  very  much  like 
Shor’s  algorithm,  in  that  it  is  a  quantum  algorithm  that  solves  the  Grover  hidden  subgroup 
problem. 


However,  this  appears  to  be  where  the  similarity  between  Grover’s  and  Shor’s 
algorithms  ends.  For  the  standard  non-abelian  QHS  algorithm  for  SN  cannot  find  the 

hidden  subgroup  Stab  ,  for  each  of  following  two  reasons: 

•  Since  the  subgroups  Stabj  are  not  normal  subgroups  of  SN ,  it  follows  from  the 

work  of  Hallgren  et  al  [55,  56]  that  the  standard  non-abelian  hidden  subgroup 
algorithm  will  find  the  largest  normal  subgroup  of  SN  lying  in  Stabj  .  But 

unfortunately,  the  largest  normal  subgroup  of  SN  lying  in  Stab,  is  the  trivial 
subgroup  of  SN . 

•  The  subgroups  Stab,,  Stab StabN_x  are  mutually  conjugate  subgroups  of S N  . 
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Moreover,  one  can  not  hope  to  use  this  QHS  approach  to  Grover’s  algorithm  to  find  a 
faster  quantum  algorithm.  For  Zalka  [77]  has  shown  that  Grover’s  algorithm  is  optimal. 


The  arguments  given  above  suggest  that  Grover’s  and  Shor’s  algorithms  are  more  closely 
related  than  one  might  at  first  expect.  Although  the  standard  non-abelian  QHS  algorithm 
on  SN  can  not  solve  the  Grover  hidden  subgroup  problem,  there  does  remain  an 
intriguing  question: 

Question.  Is  there  some  modification  of  (or  extension  of)  the  standard  QHS  algorithm 
on  SN  that  actually  solves  Grover’s  hidden  subgroup  problem? 

The  results  found  within  this  section  can  be  found  in  the  forthcoming  paper  [31]. 


II. A.  11  Distributed  quantum  computing  (DQC) 

How  can  we  use  current  or  near  future  technology  (i.e.,  technology  available  now  or 
within  the  next  five  years)  to  solve  tasks  that  we  normally  think  could  only  be  solved  by 
large  quantum  computers,  which  probably  will  not  be  available  for  at  least  twenty  to 
thirty  years  into  the  future? 


Figure  6.  Computer  architecture  model  for  quantum  distributed  computing  (QDC). 


As  an  answer  to  this  question,  we  have  in  [18,  25]  proposed  distributed  quantum 
computing  (DQC)  as  a  fast  track  roadmap  to  scalable  quantum  computing,  i.e.,  as  a 
strategy  for  effectively  using  technology  available  now  or  within  the  next  five  years  to 
perform  BIG  tasks,  normally  thought  only  possible  on  future  technology  twenty  to  thirty 
years  down  the  road. 
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II.A.12.  The  key  idea  for  DQC 

By  DQC,  we  mean  quantum  computing  on  a  network  of  small  quantum  computers 
interconnected  by  quantum  (EPR)  and  classical  channels,  as  illustrated  in  the  figure  6. 
Any  existing  or  near  future  quantum  computer  device  which,  for  example,  can  transform 
photon  (flying)  qubits  into  system  qubits  and  back  (such  as  ion  traps,  neutral  atom 
devices,  linear  optics,  etc.)  could  be  used  to  form  such  a  quantum  network. 

The  key  idea  in  [18,  25]  is  to  use  quantum  entanglement  to  distribute  control  to  the 
different  computers  within  the  above  mentioned  quantum  network.  In  particular,  our 
strategy  is  to  use  generalized  GHZ  states  to  create  cat-like  states,  which  in  turn  are  to 
be  used  to  distribute  control. 

By  a  generalized  GHZ  state,  we  mean  a  quantum  state  of  the  form 


00  0/ 


11  1, 


and  by  a  cat-like  state,  we  mean  a  quantum  state  of  the  form 


ar|00-0)  +  jff|ll-l) 

S  ' 

One  important  point  to  emphasize  is  that  once  EPR  channels  have  been  established,  each 
of  the  above  entangled  states  is  created  by  applying  only  local  unitary  operations  to  the 
individual  computing  nodes  within  the  quantum  network.  Moreover,  once  a  cat-like  state 
has  been  created,  it  does  not  matter  which  of  the  cat-like  state  qubits  is  used  for  quantum 
control.  They  all  do  operationally  the  same  thing.  This  is  illustrated  in  figure  7. 


Figure  7.  All  qubits  of  a  given  cat-like  state  are  operationally  the  same  when  it 

comes  to  quantum  control.. 
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An  n  qubit  generalized  GHZ  state  is  created  by  the  n  qubit  entangling  gate  shown  in 
figure  8, 


Figure  8.  An  n  qubit  entangling  gate. 


where  H  denotes  the  1 -qubit  Hadamard  transform,  and  where  we  have  used  the 
conventional  notation  for  the  controlled-NOT  gate. 


II.A.13.  A  universal  set  of  DQC  primitives 

Next  in  [18,  25],  we  find  a  universal  set  of  quantum  distributive  computing  primitives  for 
creating  DQC  algorithms,  and  show  systematically  how  to  convert  existing  quantum 
algorithms  into  distributed  quantum  algorithms.  (This  procedure  can  be  automated.)  In 
particular,  these  DQC  primitives  are: 

•  Cat-Creator 

•  Disentangler 

•  Reset 

•  Swap-Reset 

Cat-Creator  is  shown  in  figure  9  below, 
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Creating  a  Cat-like  State  from  a  GHZ  State 


( Cat-Creator  | 


Figure  9.  Cat-Creator  transforms  a  generalized  GHZ  state  into  a  cat-like  state. 


where  “M”  denotes  a  (standard  basis)  one  qubit  measurement  gate,  where  “X”  denotes 
the  Pauli-X  gate.  The  blue  dashed  line  indicates  that  the  classical  bit  produced  by  the 
measurement  gate  M  is  used  to  control  the  X  gates. 


The  Disentangler  primitive  is  shown  in  figure  10, 


Disentangling  the  Quantum  Channel 


«|0>|+/?|1> 


\r) 
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Second  Primitive  for  DQC 


Figure  10.  The  Disentangler  primitive  systematically  disassembles  a  cat-like  state 

after  it  has  been  utilized. 


where  “H”  again  denotes  a  one  qubit  Hadamard  gate,  “M”  a  (standard  basis)  one  qubit 
measuring  gate,  and  “Z”  a  Pauli-Z  gate.  As  indicated  by  the  blue  dashed  line,  the 
classical  bits  produced  by  the  M  gates  are  sent  to  a  classical  “logical  OR”  gate  “  ©  ” 
whose  output  is  in  turn  used  to  control  the  Pauli-Z  gate. 

The  DQC  primitives  Reset  and  Swap-Reset  are  shown  in  figures  1 1  and  12,  respectively. 
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Figure  11.  The  QDC  primitive  Reset. 


Refreshing  Channel  Qubits 


Swap-Reset 
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Fourth  Primitive  for  DQC 


Figure  12.  The  QDC  primitive  Swap-Reset. 


II.A.14.  A  systematic  procedure  for  creating  distributed  quantum  algorithms 

A  systematic  procedure  is  then  given  in  [18,25]  (which  can  be  automated)  for 
transforming  quantum  algorithms  into  distributed  quantum  algorithms.  This  systematic 
procedure  is  then  used  in  [18,25]  to  create  a  distributed  quantum  Fourier  transform,  and 
also  a  distributed  Shor  factoring  algorithm.  As  an  example,  a  distributed  quantum 
Fourier  transform  is  shown  in  figure  13. 
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Figure  13.  An  example  of  a  distributed  quantum  Fourier  transform. 


II.A.15.  The  computational  overhead  resulting  from  converting  a  quantum 
algorithm  into  a  distributed  one 

We  have  also  compared  the  complexity  of  QDC  algorithms  with  their  non-distributive 
quantum  algorithm  counterparts.  A  complexity  comparison  is  given  for  Shor’s  algorithm 
in  the  table  found  in  figure  14. 


Distributed  Factoring  Algorithm 


#  of  gates 

#  qubits 

Comm. 

Overhead 

Std  Shor 

0(«3) 

0(n) 

N/A 

Dist.  Shor 
(Theoretical) 

0(«3) 

Oin ) 

0(«3) 

Distributed 

implement 

0(n 4) 

0(n) 

Oin2) 

Figure  14.  Comparison  of  the  non-distributed  Shor  algorithm  with  the  distributed 

Shor  algorithm. 

The  asymptotic  bounds  given  in  the  above  table  are  very  conservative  estimates.  Even 
so,  it  is  surprising  and  significant  that  the  additional  time  and  space  complexity  estimates 
for  transforming  a  quantum  algorithm  into  a  distributed  one  are  essentially  washed  out  by 
the  time  and  space  complexity  of  the  shor  algorithm.  In  other  words,  for  shor’s 
algorithm,  the  additional  computational  overhead  is  insignificant!  The  same  can  be  said 
for  many  other  quantum  algorithms.  The  worst  time  complexity  overhead  we  have  seen 
resulting  from  transforming  examples  of  various  quantum  algorithms  into  distributed  ones 
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is  a  linear  slowdown  and  a  linear  increase  in  the  number  of  gates.  All  these  results 
appear  to  support  the  case  that  DQC  is  a  viable  roadmap  to  scalable  quantum  computing. 


II.A.16.  DQC  as  a  divide-and-conquer  attack  on  the  problem  of  decoherence 

We  have  also  pointed  out  that  DQC  provides  a  mechanism  for  better  dealing  with  the 
problem  of  decoherence.  It  provides  an  opportunity  for  the  application  of  a  “divide 
and  conquer”  strategy  for  dealing  with  decoherence. 

Once  EPR  channels  have  been  established,  one  need  only  focus  on  the  decoherence 
problem  for  each  spatially  separated  quantum  device  in  the  network  and  its  immediate 
environment.  The  key  idea  is  that  NOT  all  environmentally  entangling  transformations 
are  equally  likely.  In  particular,  for  spatially  separated  physical  quantum  computing 
devices,  the  most  likely  entangling  transformations  are  those  which  are  isolated  to  the 
local  quantum  device  and  its  immediate  environment.  This  is  a  substantial  simplification. 

Not  everyone  in  the  quantum  computing  community  agrees  with  this  assessment  and 
observation.  We  believe  this  stems  from  a  misunderstanding  of  what  we  have  said  above. 
What  we  are  saying  is  simply  that,  at  each  node,  one  need  only  consider  a  much  smaller 
group  (or  semigroup)  of  decohering  transformations  than  the  much  larger  group 
(semigroup)  for  the  entire  network  and  its  global  environment.  This  does  not  mean  that 
the  decoherence  at  one  node  cannot  effect  that  of  another  node,  which  is  connected  via  a 
quantum  channel. 

More  will  be  said  about  this  divide-and-conquer  approach  to  decoherence  in  future 
papers. 


II.B.  Contributions  to  the  application  of  quantum  topology  to  quantum  computing 
II.B.l.  Introduction 

In  this  section,  we  summarize  contributions  resulting  from  our  investigation  and 
exploration  of  the  application  of  quantum  topology  to  quantum  computing.  This  work 
includes  an  exploration  of  how  a  quantum  computer  could  compute  the  Jones 
polynomial,  theorems  establishing  that  generic  4x4  solutions  to  the  Yang-Baxter 
equation  are  universal  quantum  gates,  relationships  between  topological  linking  and 
quantum  entanglement,  new  universal  gates  via  solutions  to  the  Yang-Baxter  equation 
that  include  the  spectral  parameter,  new  ways  to  understand  teleportation  using  the 
categorical  formalism  of  quantum  topology,  and  a  new  theory  of  unitary  braid  group 
representations  based  on  the  bracket  model  of  the  Jones  polynomial.  These 
representations  include  the  Fibonacci  model  of  Kitaev,  and  promise  to  yield  new  insights 
into  anyonic  topological  quantum  computation. 
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Much  of  the  research  also  involves  quantum  algorithms,  and  we  have  used  functional 
integration  in  some  of  these  algorithms,  motivated  by  our  on-going  questions  about  the 
role  of  functional  integration  in  topological  quantum  field  theory.  We  are  now  working 
on  deeper  aspects  of  quantum  algorithms  associated  with  the  Jones  polynomial.  We  also 
expect  that  the  abstract  work  that  we  have  done  on  the  structure  of  teleportation  will 
impinge  on  aspects  of  distributed  quantum  computing.  Finally,  we  expect  other  insights 
from  topology  in  relation  to  the  non-abelian  hidden  subgroup  problem. 


II.B.2.  Quantum  entanglement  and  topological  entanglement 

It  is  natural  to  ask  whether  there  are  relationships  between  topological  entanglement  and 
quantum  entanglement.  Topology  studies  global  relationships  in  spaces,  and  how  one 
space  can  be  placed  within  another  (e.g.  knotting  and  linking  of  curves  in  three- 
dimensional  space).  Link  diagrams  can  be  used  as  graphical  devices  and  holders  of 
information.  In  this  vein,  Aravind  [33]  proposed  that  the  entanglement  of  a  link  should 
correspond  to  the  entanglement  of  a  quantum  state.  We  discussed  this  approach  in  [10, 
12].  Observation  at  the  link  level  is  modeled  by  cutting  one  component  of  the  link.  A  key 
example  is  the  Borommean  rings,  see  Figure  15. 


Figure  15.  Borromean  Rings. 

Cutting  any  component  of  this  link  yields  a  remaining  pair  of  unlinked  rings.  The 
Borommean  rings  are  entangled  (viz.,  the  link  is  not  split),  but  any  two  of  them  are 
unentangled.  In  this  sense,  the  Borromean  rings  are  analogous  to  the  GHZ  state 

|  GHZ)  =  (|000)  +  |lll))/V2  .  Observation  of  any  factor  (qubit)  of  the  GHZ  yields  an 

unentangled  state.  Aravind  points  out  that  this  property  is  basis  dependent,  and  we  further 
point  out  that  there  are  states  whose  entanglement  after  an  observation  is  probabilistic. 

Consider,  for  example,  the  state  (|000)  +  |00l)  +  |l0l)  +  |ll0))/2 .  Observation  in  any 

coordinate  yields  an  entangled  or  an  unentangled  state  with  equal  probability.  New  ways 
to  use  link  diagrams  must  be  invented  to  map  the  properties  of  such  states,  see  [1], 

Our  analysis  of  the  Aravind  analogy  places  it  as  an  important  question  to  which  no 
definitive  answer  has  yet  been  given.  Our  work  shows  that  the  analogy,  taken  literally, 
requires  that  a  given  quantum  state  would  have  to  be  correlated  with  a  multiplicity  of 
topological  configurations.  We  are  nevertheless  convinced  that  the  classification  of 
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quantum  states  according  to  their  correspondence  to  topological  entanglement  will  be  of 
practical  importance  to  quantum  computing,  distributed  quantum  computing  and  relations 
with  quantum  information  protocols. 


II.B.3.  Entanglement,  universality  and  unitary  R-matrices 

Another  way  to  approach  the  analysis  of  quantum  entanglement  and  topological 
entanglement  is  to  look  at  solutions  to  the  Yang-Baxter  equation  (see  below)  and 
examine  their  capacity  to  entangle  quantum  states.  A  solution  to  the  Yang-Baxter 
equation  is  a  mathematical  structure  that  lives  in  two  domains.  It  can  be  used  to  measure 
the  complexity  of  braids,  links  and  tangles,  and  it  can  (if  unitary)  be  used  as  a  gate  in  a 
quantum  computer.  We  decided  to  investigate  the  quantum  entangling  properties  of 
unitary  solutions  to  the  Yang-Baxter  equation. 

We  consider  unitary  gates  R  that  are  both  universal  for  quantum  computation  and  are  also 
solutions  to  the  condition  for  topological  braiding.  A  Yang-Baxter  operator  or  R-matrix 
[**BA**]  is  an  invertible  linear  operator  R:V  ®V  — >V  ®V ,  where  V  is  a  vector 
space,  so  that  R  satisfies  the  Yang-Baxter  equation: 

(R®  l)(l  ®  R)(R®  I)  =  (I  ®  R)(R®  l)(l  ®  R)  , 

where  /  is  the  identity  map  of  V .  This  concept  generalizes  the  permutation  of  the 
factors  (i.e.,  it  generalizes  a  swap  gate  when  V  represents  one  qubit). 

Topological  quantum  link  invariants  are  constructed  by  the  association  of  an  R  -matrix 
R  to  each  elementary  crossing  in  a  link  diagram,  so  that  an  R  -matrix  R  is  regarded  as 
representing  an  elementary  bit  of  braiding  given  by  one  string  crossing  over  another.  In 
Figure  16  below,  we  have  illustrated  the  braiding  identity  that  corresponds  to  the  Yang- 
Baxter  equation.  There  is  no  room  in  this  brief  description  to  give  the  full  translation 
from  the  topological  picture  into  the  algebraic  one.  Suffice  it  to  say  that  each  braiding 
picture  with  its  three  input  lines  (below)  and  output  lines  (above)  corresponds  to  a 
mapping  of  the  three  fold  tensor  product  of  the  vector  space  V  to  itself,  as  required  by  the 
algebraic  equation  quoted  above,  and  the  pattern  of  placement  of  the  crossings  in  the 
diagram  correspond  to  the  factors  R ®  I  and  I®R.  The  point  is  that  this  crucial 
topological  move  has  an  algebraic  expression  in  terms  of  the  R  -matrix  R  . 


Figure  16.  The  Yang-Baxter  Equation  at  the  braid  level. 
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We  worked  on  relating  topology,  quantum  computing,  and  quantum  entanglement 
through  the  use  of  R  -matrices.  In  order  to  accomplish  this  aim,  we  have  studied  the 
following  unitary  R  -matrices,  interpreting  them  as  both  braidings  and  quantum  gates. 

The  problem  of  finding  unitary  R-matrices  turns  out  to  be  surprisingly  difficult.  Dye 
(Kauffman's  former  graduate  student)  [15]  has  classified  all  such  matrices  of  size  4x4 , 
and  we  are  still  working  on  a  general  theory  for  the  classification  of  unitary  R  -matrices 
in  other  dimensions. 

A  key  question  about  unitary  R-matrices  is  to  understand  their  capability  of  entangling 
quantum  states.  We  use  the  criterion  that  |  (p)  -  a  1 00}  +  b |  Ol)  +  c 1 10)  +  d 1 1 1)  is  entangled 
if  and  only  if  ad  -be  ^0  .  This  criterion  is  generalized  to  higher  dimensional  pure  states 
in  the  papers  [10,  11,  12,  20]  by  Kauffman  and  Lomonaco.  We  discovered  families  of  R  - 
matrices  that  detect  topological  linking  if  and  only  if  they  can  entangle  quantum  states.  A 
recent  example  in  [16]  is  a  unitary  R  -matrix  that  is  highly  entangling  for  quantum  states. 
It  takes  the  standard  basis  for  the  tensor  product  of  two  single-qubit  spaces  onto  the  Bell 
basis.  On  the  topological  side,  R  generates  a  non-trivial  invariant  of  knots  and  links  that 
is  a  specialization  of  the  well-known  link  invariant,  the  Homflypt  polynomial. 

Entanglement  and  quantum  computing  are  related  in  a  myriad  of  ways,  not  the  least  of 
which  is  the  fact  that  one  can  replace  the  CNOT  gate  by  another  gate  R  and  maintain 
universality  (as  described  above)  just  so  long  as  R  can  entangle  quantum  states.  That  is, 
R  can  be  applied  to  some  unentangled  state  to  produce  an  entangled  state.  It  is  of  interest 
to  examine  other  sets  of  universal  primitives  that  are  obtained  by  replacing  CNOT  by 
such  an  R . 

We  proved  that  certain  solutions  R  to  the  Yang-Baxter  equation,  together  with  local 
unitary  two  dimensional  operators,  form  a  universal  set  of  quantum  gates.  Results  of  this 
kind  follow  from  general  results  of  the  Brylinskis  [39]  about  universal  quantum  gates. 

The  Brylinskis  show  that  a  gate  R  is  universal  in  this  sense,  if  and  only  if  it  can  entangle 
a  state  that  is  initially  unentangled.  We  show  that  generically,  the  4x4  solutions  to  the 
Yang-Baxter  equation  are  universal  quantum  gates. 

For  example,  the  following  solutions  to  the  Yang-Baxter  equation  are  universal  quantum 
gates  (in  the  presence  of  local  unitary  transformations): 
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where  a,  b  are  unit  complex  numbers  with  a2  ^b2 . 
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The  unitary  matrix  if  is  the  Bell-Basis  change  matrix,  alluded  to  above.  The  unitary 
matrix/?'  is  a  close  relative  to  the  swap-gate  (which  is  not  universal).  The  unitary 
matrix/?"  is  both  a  universal  gate  and  a  useful  matrix  for  topological  purposes  (it  detects 
linking  numbers).  In  this  last  example,  we  see  a  solution  to  the  Yang-Baxter  equation 
that  detects  topological  linking  exactly  when  it  entangles  quantum  states. 

These  results  about  R  -matrices  are  fundamental  for  understanding  topological 
relationships  with  quantum  computing,  but  they  are  only  a  first  step  in  the  direction  of 
topological  quantum  computing.  In  topological  quantum  computing  one  wants  to  have 
all  gates  and  compositions  of  gates  interpreted  as  part  of  a  single  representation  of  the 
Artin  Braid  Group.  By  taking  only  a  topological  operator  as  a  replacement  for  CNOT , 
we  leave  open  the  question  of  the  topological  interpretation  of  local  unitary  operators. 

One  must  go  on  and  examine  braiding  at  the  level  of  local  unitary  transformations  and  the 
problem  of  making  fully  topological  models.  The  first  step  in  this  process  (although 
made  only  recently  by  us  [30])  is  to  classify  representations  of  the  three-strand  braid 
group  into  5/7(2) .  Here  one  looks  not  for  solutions  to  the  Yang-Baxter  equation,  but 
rather  for  matrices  A  and  B  in  U (2)  such  that  ABA  =  BAB .  This  is  the  analogue  of 
the  Yang-Baxter  equation  and  it  has  many  solutions  in  Z7(2)  .  Some  of  these  pairs 
generate  dense  subsets  of  U (2)  and  so  can  be  used  in  principle  to  generate  all  local 
unitary  transformations.  At  this  stage  one  has  a  mixed  topological  generation  of  quantum 
computing:  one  type  of  transformation  for  local  unitary  operators  and  a  second  (Yang- 
Baxter)  transformation  for  the  universal  2-qubit  gate. 


evolution 


Figure  17.  A  topological  quantum  computer. 


In  the  diagram  above,  we  have  illustrated  the  form  of  such  a  rudimentary  topological 
quantum  computer.  The  crossings  represent  2-qubit  gates  that  are  solutions  to  the  Yang- 


38 


Baxter  equation.  The  dark  circles  represent  local  unitary  transformations  that  may 
themselves  be  generated  by  a  unitary  representation  of  the  braid  group.  To  go  further 
involves  finding  braiding  representations  into  U (2)  that  extend  to  dense  representations 
in  U(N)  for  larger  values  of  N  .  This  is  where  topological  quantum  field  theory  comes 
into  play. 

In  the  next  section  we  outline  our  approach  to  full  topological  quantum  computation. 


II.B.4.  Topological  quantum  Held  theory  and  topological  quantum  computation 

As  described  above,  one  comes  to  a  barrier  if  one  only  attempts  to  construct  individual 
topological  gates  for  quantum  computing.  In  order  to  go  further,  one  must  find  ways  to 
make  global  unitary  representations  of  the  Artin  Braid  Group.  One  way  to  accomplish 
this  aim  is  via  topological  quantum  field  theory.  Topological  quantum  field  theory 
originated  in  the  work  of  Witten  [73]  with  important  input  from  Atiyah  [34],  This  work 
opened  up  quantum  field  theoretic  intepretations  of  the  Jones  polynomial  (an  invariant  on 
knots  and  links,  new  at  that  time)  and  gave  rise  to  new  representations  of  the  braid 
groups.  The  basic  ideas  of  topological  quantum  field  theory  generalize  concepts  of 
angular  momentum  recombination  in  classical  quantum  physics.  In  [65,  30],  we  use 
generalizations  (so-called  q-deformations)  of  the  Penrose  [68]  formalism  of  spin 
networks  to  make  models  of  topological  quantum  field  theories  that  are  finite 
dimensional,  unitary  and  that  produce  dense  representations  of  the  braid  group  into  the 
unitary  group.  These  representations  can  be  used  to  do  quantum  computing.  In  this  way, 
we  recover  a  version  of  the  results  of  Freedman  [45,  46,  47,  48,  49]  and  his  collaborators 
and,  by  making  very  concrete  representations,  open  the  way  for  many  applications  of 
these  ideas.  Our  methods  are  part  of  the  approach  to  Witten's  invariants  that  is 
constructed  in  the  book  of  Kauffman  and  Lins  [65].  This  work  is  directly  based  on  the 
combinatorial  knot  theory  associated  with  the  Jones  polynomial.  Thus  our  work  provides 
a  direct  and  fundamental  relationship  between  quantum  computing  and  the  Jones 
polynomial. 

Here  is  a  very  condensed  presentation  of  how  unitary  representations  of  the  braid  group 
are  constructed  via  topological  quantum  field  theoretic  methods.  The  structure  described 
here  is  sometimes  called  the  Fibonacci  model  [30,  60,  69],  One  has  a  mathematical 
particle  with  label  P  that  can  interact  with  itself  to  produce  either  itself  labeled  P  or 
itself  with  the  null  label  * .  When  *  interacts  with  P  the  result  is  always  P.  When  * 
interacts  with  *  the  result  is  always  * .  One  considers  process  spaces  where  a  row  of 
particles  labeled  P  can  successively  interact  subject  to  the  restriction  that  the  end  result  is 
P . 

For  example  the  space  F[(<iA)cJ  denotes  the  space  of  interactions  of  three  particles 
labeled  P  .  The  particles  are  placed  in  the  positions  a,b,c .  Thus,  we  begin  with  (PP)P . 
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In  a  typical  sequence  of  interactions,  the  first  two  Ps  interact  to  produce  a  * ,  and  the  * 
interacts  with  P  to  produce  P  . 

(PP)P  ->(*)P  ->  P . 

In  another  possibility,  the  first  two  P  ’s  interact  to  produce  a  P ,  and  the  P  interacts 
with  P  to  produce  P  . 

(PP)P  -»(P)P  ->  P . 

It  follows  from  this  analysis  that  the  space  of  linear  combinations  of  processes 
F[(«6)c]  is  two  dimensional.  The  two  processes  we  have  just  described  can  be  taken 
to  be  the  qubit  basis  for  this  space.  One  obtains  a  representation  of  the  three  strand  Artin 
braid  group  on  F[(aZ>)c]  by  assigning  appropriate  phase  changes  to  each  of  the 

generating  processes.  One  can  think  of  these  phases  as  corresponding  to  the  interchange 
of  the  particles  labeled  a  and  b  in  the  association  (ab)c .  The  other  operator  for  this 
representation  corresponds  to  the  interchange  of  b  and  c . 


This  interchange  is  accomplished  by  a  unitary  change  of  basis  mapping 

F  :V[(ab)c]^v[a{bc)]. 

If 

A :  F[(aft)c]  — >  F[a(6c)] 

is  the  first  braiding  operator  (corresponding  to  an  interchange  of  the  first  two  particles  in 
the  association),  then  the  second  operator 

B  :F[(«&)c]->f[«(Ac)] 


is  accomplished  via  the  formula  B  =  F " 1 AF  where  the  A  in  this  formula  acts  in  the  second 
vector  space  F[(a6)c]  to  apply  the  phases  for  the  interchange  of  b  and  c . 


In  this  scheme,  vector  spaces  corresponding  to  associated  strings  of  particle  interactions 
are  interrelated  by  recoupling  transformations  that  generalize  the  mapping  F  indicated 
above.  A  full  representation  of  the  Artin  braid  group  on  each  space  is  defined  in  terms  of 
the  local  interchange  phase  gates  and  the  recoupling  transformations.  These  gates  and 
transformations  have  to  satisfy  a  number  of  identities  in  order  to  produce  a  well-defined 
representation  of  the  braid  group.  These  identities  were  discovered  originally  in  relation 
to  topological  quantum  field  theory.  In  our  approach  [30],  the  structure  of  phase  gates  and 
recoupling  transformations  arise  naturally  from  the  structure  of  the  bracket  model  for  the 
Jones  polynomial.  Thus  we  obtain  a  knot-theoretic  basis  for  topological  quantum 
computing. 

Many  questions  arise  from  this  approach  to  quantum  computing.  The  deepest  question  is 
whether  there  are  physical  realizations  for  the  mathematical  particle  interactions  that 
constitute  such  models.  It  is  possible  that  such  realizations  may  come  about  by  way  of 
the  fractional  quantum  Hall  effect  or  by  other  means.  We  are  working  on  the  physical 
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basis  for  such  models  by  addressing  the  problem  of  finding  a  global  Hamiltonian  for 
them,  in  analogy  to  the  local  Hamiltonians  that  can  be  constructed  for  solutions  to  the 
Yang-Baxter  equation.  We  are  also  investigating  specific  ways  to  create  and 
approximate  gates  in  these  models,  and  we  are  working  on  the  form  of  quantum 
computers  based  on  recoupling  and  braiding  transformations. 

These  models  are  based  on  the  structure  of  the  Jones  polynomial  [4,  61,  62,  63,  64], 

They  lead  naturally  to  the  question  of  whether  or  not  there  exists  a  polynomial  time 
quantum  algorithm  for  computing  the  Jones  polynomial.  The  problem  of  computing  the 
Jones  polynomial  is  known  to  be  classically  #P-hard,  and  hence,  classically 
computationally  harder  than  NP-complete  problems.  Should  such  a  polynomial  time 
quantum  algorithm  exist,  then  it  would  be  possible  to  create  polynomial  time  quantum 
algorithms  for  any  NP-complete  problem,  such  as  for  example,  the  traveling  salesman 
problem.  This  would  indeed  be  a  major  breakthrough  of  greater  magnitude  than  that 
arising  from  Shor's  and  Simon's  quantum  algorithms.  The  problem  of  determining  the 
quantum  computational  hardness  of  the  Jones  polynomial  would  indeed  shed  some  light 
on  the  very  fundamental  limits  of  quantum  computation. 

A  polynomial  time  quantum  algorithm  (called  the  Ahronov-Freedman-Jones-Kitaev- 
Landau  (AFJKL)  algorithm)  for  approximating  the  value  of  the  Jones  polynomial  L(t) 
at  primitive  roots  of  unity  can  be  found  in  [32],  We  are  currently  writing  a  paper  [30]  that 
shows  that  this  algorithm  can  not  successfully  be  extended  by  polynomial  interpolation  to 
a  polynomial  time  quantum  algorithm  for  computing  the  Jones  polynomial.  However, 
there  is  a  loop  hole.  It  may  well  still  be  possible  to  modify  the  AFJKL  algorithm  in  such 
a  way  that  it  can  be  used  to  create  a  polynomial  time  algorithm  for  L(t) .  We  propose  to 
investigate  why  this  is  or  is  not  the  case.  Our  objective  is  to  come  to  a  better 
understanding  of  the  exact  divide  between  classical  and  quantum  algorithms. 


III.  Epilogue 

In  the  above  section  II,  we  have  explained  extensively  the  many  contributions  to  quantum 
computing  and  to  quantum  information  science  resulting  from  this  contract,  as  they  were 
outlined  in  the  executive  summary  given  in  section  I.  Supporting  documents  and 
supplementary  material  for  this  final  report  can  be  found  at  the  website: 

http://www.csee.umbc.edu/~lomonaco/DARPA/01-06fmalrpt 
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